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تحلیل محبط نیمه بی‌نهایت همسان جانبی به پیچش وارد بر جداره‌ی حفره با طول محدود* 


محمدرضا ۱۳ مرتضی اسکندری قادی 

چکیده در این مفاله یک محیط نیمه بی‌نهایت با رفتار ایزوتروپ (همسان) جانبی که محور ایزوتروپی (همسانی) آن عمود بر سط حآزاد 
است و حفره‌ی استوانه‌ای با طول محدود در امتداد محور ایزوتروپی د رآن ایجاد شده است, در نظ رگرفته شده و پاس خحآن به پیچش معلوم 
روی دیواره‌ی این استوانه و حول محور استوانه به صورت تحلیلی بررسی می‌شود. بدین منظور معادلات تعادل استاتیکی حاکم بر مساله در 
دستگاه محتصات استوانه‌ای نوشته شده و با تقسیم محیط به دو ناحیه و نوشتتن معادلات تعادل برای هر ناحیه به صورت مجزا و استفاده از 
تبدیل کسینوسی فوریه, جابه‌جایی محیط در فضای تبدیل‌یافته ارائه می‌گردد. به کمک قضیه عکس تبدیل انتگرالی جابه‌جایی‌ها در فضای 
اصلی مسنله حاصل می‌شود. با نوشتن شرایط مرزی و پیوستگی, معادله‌ی انتگرالی کوشی حاکم بر مسأله به دست می‌آید. با حل معادله‌ی 
انتگرالی حاکم, توابع تنش و تغییر مکان در هر نقطه از محیط به دست می/یند. نتایج به دست آمده برای محیط‌های همسان جانبی با نتایج 


موجود برای محیط‌های همسان مفایسه می‌شود. 


واژه‌های کلیدی پیچش. همسان جانبی نیم‌فضاء حفره‌ی استوانه‌ای» تبدیل کسینوسی فوریه معادلات انتگرالی. 
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* تاریخ تصویب مقاله ۸٩/۱۱/۰۲‏ و تاریخ دریافت نسخه نهایی اصلاح شده ٩۰/۸/۱۰‏ 


)۱( دانشجوی دکترا مهندسی عمران- سازه. گروه مهندسی عمران. دانشکده فنی و مهندسی. دانشگاه علم و فرهنگ. 
(۲) نویسنده‌ی مسوژول: دانشیار دانشکده مهندسی عمران. پردیس دانشکده‌های فنی» دانشگاه تهران. 
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مسأله‌ی انتشار امواج در یک محیط ناشی از بارگذاری 
خارجی و تحلیل پاسخ یک محیط الاستیک به اعمال 
تنش بر روی جدار حفره‌ای در داخل آن از جمله 
مباحثی است که در قرن گذشته توجه بسیاری از 
محققان و مهندسان ریاضیات کاربردی و مکانیک 
مهندسی را به خود جلب کرده است. مسأله‌ی بررسی 
میزان تتش و کرنش به وجود آمده در یک محیط 
نیمه‌بی‌نهایت تحت اثر نیروهای ناشی از شمع‌های 
کوتاه يا بلند نیازمند به بررسی اندرکنش شمع و محیط 
است و این خود نیاز به تحلیل محیط نیمه بی‌نهایت 
تحت اثر نیروهای ناشی از شمع موثر بر فصل مشترک 
جدار شمع و محیط دارد. 

وسترگارد. اولین بار در سال ۰۱۹۶۱ یک محیط 
بی‌نهایت با رفتار ارتجاعی با یک حفره‌ی استوانه‌ای با 
طول بی‌نهایت را در نظر گرفت و پاسخ استاتیکی آن به 
فشار شعاعی روی طول محدودی از دیواره‌ی حفره را 
به دست آورد [15]. آنالیز سه‌بعدی حفره‌ی استوانه‌ای 
توسط گرین و زرنا در سال ۱۹6۶ انجام گرفت [3] در 
تل ۱۹۵۲ سیر که به تعلیل بل شخیط ارتطاقی 
تحت اثر فشار آنی بر سطح داخلی حفره پرداحت 
[12] بعدها حردن دی سال ۱۹۲۲ تال مشاه .ربا 
اعمال فشار دینامیکی خطی بر روی قسمتی از دیواره‌ی 
حفره حل کرد [4]. در سال ۱۹۸۰ پارنز به مطالعه‌ی 
اعمال تنش پیچشی خطی متحرک که به صورت 
دایره‌ای در امتداد محور داخلی حفره‌ی استوانه‌ای» واقع 
در یک محیط الاستیک بی‌نهایت در حرکت بود. 
پردانعت [10]. هم‌چنین اسنیدون در سال ۱۹۵۲ و کل 
و هاث در سال ۱۹۵۸ حالت دو بعدی یله پارنز 
را مورد بررسی قرار داده بودند [1,14]. پارنز در سال 
ویس اهوال سا ماع وی عفن 
خطی, به طور مجزاء. بر روی دیواره‌ی حفره‌ی نامحدود 
در یک محیط بی‌نهایت الاستیک پرداخت و نتایج 
عددی برآورد ملفه‌های تنش و جابه‌جایی در نزدیکی 
نیرو و در طول حفره را ارائه نمود [7]. پارنز در سال 


تشریه مهندسی عمران دانشگاه فردوسی مشهد 


تحلیل محیط نیمه بی نهایت همسان جانبی به پیچش وارد بر.. 


۳ مسأله‌ی پیچش خود را در حالت دینامیکی حسل 
کرد و نتایج آن را با حالت استاتیکی قبلی؛ مقایسه 
نمود [8]. در اين مقاله جابه‌جایی‌ها و تنش‌های حاصل 
از اعمال فشار خطی پر جداره‌ی حفره‌ی استوانه‌ای به 
صورت انتگرال‌های بی‌نهایت بیان شده‌اند. 

با توجه به بررسی‌های انجام شده در نزدیکی 
عمال نیروی متمرکز, همگرایی انتگرال‌ها پسه سرعت 
کاهش می‌یابد و نقاط تکین خود را نشان می‌دهند. از 
ین‌ری نتایج در نقاطی دور از نقاط تکین دارای دقت 
خوبی می‌باشند. این محقق در سال ۱۹۸1 پاسخ محیط 


حفره‌ی استوانه‌ای با طول بی‌نهایت را به دست آورد 
[11]. وجود حفره با طول محدود در محیط‌های 
نیمه‌بی‌نهایت. به علت وجود سختی محیط در کف 
حفره. باعث دشواری تحلیل مسأله هی گردده از اینن‌ره 
در اکثر تحقیقات اثر اعمال نیرو بر جداره‌ی حفره‌ی 
نامحدود در یک محیط مورد مطالعه قرار گرفته است. 
از آن‌جا که در بسیاری از کاربردهای مهندسی از جمله 
حفر فونداسیون‌هاء مطالعات ژئوتکنیکی در محل و 
مدل کردن گسترش تنش وارد بر دیواره‌ی حفره با 
حفره‌ای با طول محدود سر و کار داریم؛ مطالعه‌ی اثر 
نیرو بر حفره با طول محدود مهم می‌باشد. 

پاک و انا آوم وی بان ۱۹4۲ ان کنشن هی شر 
دیواره‌ی حفره‌ای محدود در یک محیط همگن و 
همسان را با دقت مورد بررسی قرار داده‌اند [6]. در 
حال حاضر با توجه به استفاده‌ی روزافزون از مواد 
ناهمسان نیاز به تحقیقات در زمینه‌ی انتشار امواج در 
این محیط‌ها بیشتر احساس می‌شود. برای مثال. مواد 
کامپوزیت که در سال‌های اخیر در زمینه‌ی علوم 
میس کارت سای هت دای هی ی 
ناهمسانی می‌باشند. از سوی دیگر در زمین‌هایی که 
خاک تحت اثر نیروی ثقلی رسوب کرده است و 
نهشته‌های طبیعی سربارشده روی هم تشکیل داده 
است. خاصیت ناهمسانی وجود دارد؛ اما با توجه به 
ملاحظات کاربردی در زمینه‌ی مهندسی. محیط‌های 
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ناهمسان معمولاً به صورت همسان جانبی و با 
ارتوتروپیک مدل‌سازی می‌شوند [۱71]. در حالت کلی 
در صورتی که در یک محیط سه‌بعدی. فرض همسان 
بودن ماده کنار گذاشته شود تحلیل تیش از پیچیده و 
طولانی خواهد شد. 

در این مقاله یک محیط نیمه بی‌نهایت با رفتار 
همسان جانبی که حفره‌ی استوانه‌ای با طول محدود در 
امتداد محور ایزوتروپی در آن ایجاد شده, در نظر گرفته 
شده است و پاسخ آن به پیچش معلوم روی دیواره‌ی 
این استوانه و حول محور استوانه به صورت تحلیلی 
بررسی می‌شود. بدین منظور معادلات حرکت حاکم بر 
مسأله در دستگاه مختصات استوانه‌ای نوشته می‌شود. با 
تقسیم محیط به دو ناحیه و نوشتن معادلات حرکت 
برای هر ناحیه به صورت مجزا و برقراری شرایط 
پیوستگی و مرزی پاسخ محیط تعیین می‌شود. در روش 
حل از تبدیل کسینوسی فوریه استفاده شده است؛ نذا 
جابه‌جایی‌های محیط در فضای تبدیل‌یافته به دست 
هی آبشتن: اه کشک قبت هی عکسی تسیل انیک اس + 
جابه‌جایی‌ها در فضای اصلی مسأله حاصل می‌شود. 
ثابت‌های انتگرال‌گیری با معرفی یک تابع میانی و 
ارضای شرایط پیوستگی و سازگاری بین دو ناحیه در 
فضای واقعی به دست مس آیتتد: با نوشتن شرایط 
سازگاری. معادله انتگرالی کوشی حاکم بر تابع میانی 
مسأله حاصل می‌شود. با حل معادله‌ی انتگرالی کوشی, 
جابه‌جایی‌ها به دست می آید و با استفاده از معادلات 
رفتاری تنش‌ها در محیط همسان جانبی به دست 
شرا با برآورد عددی نشان داده می‌شود که 
جواب‌های مسأله در حالت ساده‌تن مربوط به 
محیط‌های همسان. بر جواب‌های موجود منطبق است 


که دلالت بر صحت و دقت نتایج دارد. 


معادلات حاکم بر مسأله مقدار مرزی 


یک محیط نیمه بی‌نهایت ارتجاعی همگن و با رفتار 
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۳۷ 


ای ی ای تایب یت 1 
شعاع 2 مطابق شکل (۱) در نظر گرفته می‌شود. 
دستگاه مختصات استوانه‌ای (0,2,) را مطابق این 
شکل در نظر می‌گيریم. محور ایزوتروپی محیط در 
امتداد محور 7 و در نتيجه عمود بر سطح نیم فضا در 
نظر گرفته می‌شود. دیواره حفره در فاصله‌ی 1> 7 > 0 
ف اردشت پرشی مر (ع راز درد میتی 
افش تین یی کستعاور معلوم هرد 
می‌باشد. به منظور حل مسأله» محیط را به دو ناحیه‌ی 


یک و و به صورت ذیل تقسیم می‌کنیم (شکل 0 
)۱( 01 < 2ر22 > 0 > 0,ه < ۲ |(6,0,2 ) < بظ 


)۳ (1< 22,2 >0 > 0,ه > ۲ |(60,2) )< با 


2) 
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ی 
شکل ۱ حفره‌ی استوانه‌ای در محیط نیمه بی‌نهایت 


با توجه به این تقسیم‌بندی, شرایط مرزی مسأله به 


,۳( 0621 ,() ۲ (2,2)ر 
3 ۲<6 ,0 (۲,۵) رو 
62 > ,0 2 (1و) رورت 
)1( ۲<3 0-0 2 

02 
ِ‌ > ۹2 


,20 (۲,1) 2 
‌ مج 
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و شرایط در بی‌نهایت به صورت زیر نوشته می‌شود: 
مب ( 2+ کل 


در این روابط اندیس ۱ مربوط به ناحیه‌ی ,5 و 


)۸ 2ب [ 0 (7,) بل 


اندیس ۲ مربوط به ناحیه ی می‌باشد. به علاوه در 
این روابط وت و ور تنش‌های برشی و تنها تتش‌های 
غیر صفر محیط هستند و هم‌چنین ولا < لا تغییرمکان 
در امتداد 0 در دستگاه مختصات معرفی شده است. 
طبیعت مسأله چنان است که ۶0 رل رلا منی‌باشند. 
روابط بین تنش - تغییرمکان برای تنش‌های غیر صفر 
به صورت زیر می‌باشند [د]: 

1 


)( ۹) 


5 «لا 1 
9۹ 1 05 ۳ 


افو 
99 )۳ و 


که در آن ۲ مدول برشی در صفحه ایزوترویی و ۳ 
مدول برشی در صفحه عمود بر صفحه ایزوترویی 
می‌باشد. معادله‌ی تعادل بر حسب تغییرمکان در هر 

67 ء بل _ یلا۵ 1 02 


2۳۵-2 9 22 5 ۱۱) 


5 ض‌ ر لا 1202 تا 0 
22 ۴ ۲ ۲ ۵۸8 


که در آن ی ۳ سبت سختر پرشی در صفحه‌ی 
۳ 

ایزوتروپی به سختی برشی در صفحه‌ی عمود بر 

ضس‌خاای: ایزو ترویسی نی باشسانه. پیوسستکی نش ق 

تغییرمکان در فصل مشترک ناحیه‌های ,5 و رل به 


(۱۳( 1< 2 و( 2 29) و۲ < (2,2) رم 
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ره (ه رم 2۳ 
,7 


<1 ۱۶( 
94 ۲: ۱ 


(2,2) رور؟ < (2,2) ,ور برای 1< 2 است. 


حل معادلات حاکم 
با توجه به دامنه‌ی مسأله برای تبدیل معادلات با 
مشتقات جزئی (۱۲ و ۱۱) به معادلات دیفرانسیل 
معمولی از تبدیل کسینوسی فوریه استفاده می‌کنیم 
[13]. از آن‌جایی که دامنه‌ی ,8 از 2-0 تا بی‌نهایت 
و دامنه‌ی ی از 2<1 تابی‌نهایت است. تبدیل 
کسینوسی فوریه توابع در اين دامنه‌ها متفاوت است و 


به ترتیب به صورت زیر نوشته می‌شوند: 


(۵ 7 (2ع )09 (2 "۳ ]2 (ع ,0 
۳۹ مه( ماهر (٩,‏ 


که در آن ,؟ و ؟ توابع دل‌خواه به ترتیب با دامنه‌های 
5 و ی هستند و بل و تبدیل فوریه کسینوسی 
آنها می‌باشند. بنابراین» تبدیل کسینوسی فوریه معادلات 
(۱۲ و ۱۱) به صورت زیر در می‌آید: 

1 اه 


1 
سس 
1 


۱۳۳ 11 
(۱۷) 2را< 1 0- بت( + تا 


که.در آن از شرایط (1۳۸) اسفاده شبده اسشت. با تقییر 
متغیر یا <ی1 و 06-6 معادله‌ی (۱۷) به معادله‌ی 
بسل اصلاح‌شده‌ی مرتبه‌ی اول در می‌آید و جواب 
آن‌ها بر حسب 7 و چّ در حالت کلی به صورت: 


(021) ,کا(ج) ,۸ ع (و,۲) رل 


(۸ 1,2 - زٌ ,(05۲) ,آ(ج) ,۲ + 


سال بیست و چهارم. شماره یک, ۱۳۹۱ 
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اصلاح شده‌ی نوع اول و دوم از مرنبه‌ ی یک می‌باشند. 
در ناحیه‌ی ,۲ هنگامی که 0ب ۲ رآ به سمت 
قبول (8),ظ صفر می‌باشد. هم‌چنین در ناحیه‌ی 8 
هنگامی که ۲-0 ,16 به بی‌نهایت میل می‌کند» در 
این صورت () بش۸ صفر می‌باشد. با استفاده از 
قضیه‌ی عکس تبدیل کسینوسی فوریه جابه‌جایی‌ها در 


دو ناحیه به ترتیب به صورت: 


رکب ]هب 


(۱۹) 
*005)( 1 ۲ <8,2,<0 

و ریم باقابظ ‏ | - ۵ب 
1< 2 و > ۲ 5 (1- 605602« 


تلا دستست فیس تال تیسا دون نطتر ترفن شبرایط 
پیوستگی(۱۳) روی مرز استوانه‌ای شکل در ناحیه ,3 


می توان ,مت را به صورت: 


۲7 0( >> 
10( 221 


01 
)۳( (ه<() مرچ 


و تابع ووت در مرز ناحیه بارا به شکل: 


(۲۲( 1< 7 ,(۱2/)2 - (2 )تس هن 
معرف تابع کرنش پرشی در صفحات موازی صفحه‌ی 
۲-0 در ۲<8 بر حسب عمق می‌باشد. با جاگذاری 
روابط (۲۰ و )۱٩‏ در (۲۲ و ۱) و استفاده از روابط 
مشتق توابع بسل اصلاح شده و هم‌چنین عکس تبدیل 
کسینوسی فوریه ضرایب (۸,)5 و ()3 به صورت 
ذیل به دست مش اقلا 


(0۲۳ 427 (۲)2(۰05)62 ۳ < (و) ,۸ 


سال بیست و چهارم شماره یک ۱۳٩۱‏ 


۳۹ 


ع عم (ه)ر ی - (8),ظ 
با جاگذاری روابط (۲۶ و ۲۳) به ترتیب در روابط 
(۲۰ و )۱٩‏ میزان جابه‌جایی در دو ناحیه بر حسب 
تابع مجهول (7/2 به دست می‌آید. جابه‌جایی‌های نا 
و ولا باید شرط پیوستگی (۱۶) را ارضا نمایند. با 
مشتق‌گیری از ,لا و ولا نسبت به 2 داریم: 
0,۳۶9 ال - رم 


(۲۵) ((و + 2 ,۲),ه + 
0 < 7 وه < ۲ ,ی (ی) < 


۳ ۳۳ ۳ 
 ):-‏ اخ - (هر خ2 


((21 - 4 + ۲,2),ه + 


1< 72,ه > ۲ ری (ی) « 
)۳۲۹( 
به طوری که 
‌ ۳ - ۵,0,۵ 
0 0 2 5 0,06 


و با جاگذاری (۲۳ و ۲۵) در شرط پیوستگی (۱۶» این 
شرط به صورت معادله انتگرالی زیر برای تعیین (5)2 


با 1< 72 به دست می‌آید: 
2-20 ۵062-۵۵6 | 
هس66 ]+ 


اجه یق(0با+ست 6 ]- 


معادله‌ی (۲۹) معادله انتگرالی تعمیم‌یافته‌ی کوشی 


نامیده می‌شود. 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


حل معادله‌ی انتگرالی حاکم 
همان‌طور که در معادله‌ی (۲۹) مشاهده می‌شود توابع 
و ,ه باید برای 8<- ۲ بررسی شوند. با محاسبه‌ی 
(۵, "۵,62 و (۵, 4,)2 و تعریف متغیرهای بی‌بعد 
_ 2 


0 ٍ 2 7 معادله‌ی (۲۹) به صورت 
2 2 2 


واه هی شنود: 
مر مرح 1 1 2« 
6 توا | 
2ب 2-2ببا سپ 1 


مه( ]+ 


2<1 (80-:2(004+ )6 - 
به طوری که 
(۳۱( (702 - (20 
ی 
۳ و | - (۵,! 
ی ایک _ 
۳۳ ا وج ,| ۵ 


(مه- مقر خه)رهاه ] - 10 
4 (ه) > [(ه + 2۵,(ه),ه + 


(۳( 


با حل معادله‌ی انتگرالی (۳۰) تنش‌های برشی 
کشرز ۱8 ان لح هه اذمست:می ینت سرا ضرع 
معادله‌ی انتگرالی (۳۰) از تغییر متغیرهای زیر به منظور 
بی‌بعد کردن استفاده می‌کنیم: 


1/6 ۲ ۰1/2 )۳۵( 


به کمک رابطه‌ی  )۳۵(‏ معادله‌ی (۳۰) را می‌توان 


9 


تشریه مهندسی عمران دانشگاه فردوسی مشهد 


تحلیل محیط نیمه بی نهایت همسان جانبی به پیچش وارد بر.. 


موه و6002 [ 
0 


نوشت» به طوری که: 


ای ۰ 

۳ ی 7 7 
[رب ۱ ٩‏ ,+( 1 0 
۷ ۷ ۷ 

5 

2)( - ۶0 /۰( (۳ 


با توسیع تابع 20/0 به صورت زوج نسبت به 
مبداً می‌توان رابطه‌ی (۳۳) را به صورت زیر نوشت: 


0 ۳ 


۷ 


‌ 


( 201 /۳( ۵۷ - 


با بررسی رفتار تابع (262 در نقاط تکینء می‌توان 
تابع را در فاصله‌ی [ابل-] به صورت: 


(6۰ ترفن ومع(« 20 


توت که دز آن تابع مجهول (10 یک تابع هموار در 


(۰ انتگرال (۳۹) به صورت زیر بیان می‌شود: 


۷ و2۶ 


‌ 


0 


1 1 
2 ۷ - ( 


0-۷ ۱( 


با در نظر گرفتن نقاط تکین انتگرال فوق در 
محدوده‌ی [لبل-] و استفاده از بسط به صورت حند 
جمله‌ای ژاکوبی (() /!۳) با << برای تابع 
تکین موجود در انتگرال و نیز زوج بودن تابع 100 
برای 1> ۷ > 1-. معادله‌ی انتگرالی به صورت دستگاه 


معادلات زیر نوشته می‌شود [2: 
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/ ۳ 


5 


1 21 
2 


2-1رساحز ‏ وار)8 ,۲6۳ 


که در آن از روش تجمع محلی (۱166۵0 عمنایم0011) 
[۱۷] برای نقاط تکین تابع زير علامت انتگرال در (4۱) 
استفاده شده است. در رابطه‌ی (۶۲): 


(بمب+م+۵ +( 
(۳: +۵ +2 1(۲+/۲12(۲+۵ 
۳( جم +م +1021 


۹ 


و ,۷,۶ ريشه های چند جمله‌ای ژاکوبی 


ص 

(61) 2۲-1 ...با ز ,0 (۳) ورط 
رس 

(1۵) اه ا یا بط 


می‌باشند. از آن‌جایی که رل و ,۷ نسبت به مبداً 
متقارن هستند. دستگاه معادلات (4۲) برای [....,1-ز 


باه بو ریت 
آ 
ار 
1-1 


و برای 20 ,, به صورت زیر نوشته می‌شود: 


آ 
)۷ 0 2۷ 


رابطه‌ی (47) شامل 3 معادله می‌باشد که با حل 
هم‌زمان آن‌ها (,10۷ برای "...1,۰ به دست می‌آید. 
با توجه به رابطه بین (16 و (5)2 و هم‌چنین روابط 
(۲۶ و ۲۳ و ۲۰ و )۱٩‏ پاسخ محیط به تنش پیچشی 
معلوم "7 قابل محاسبه می‌باشد. در بسیاری از 


کاربردهای مهندسی تنش و جابه‌جایی روی دیواره‌ی 


سال بیست و چهارم شماره یک ۱۳٩۱‏ 


3 


حفره مورد نیاز است. از این‌ری با استفاده از رابطه 


() می‌توان جابه‌جایی را در 8< ۲ به صورت زیر 
محاسبه نمود: 


(002)- دک - (2 ,)۱ 
(۸ ۱ ۱ ۱ ۱ 0 ۲ 
(46 1( - 8 ,و +( +9,)2] - | + 


به طوری که: 


اکز ]_ _ 
٩‏ وی وج , [-- 9:6۵ 


دم 0ج 6-0,و++یما ]-96 
(۵۱) (0 ۲ - (6 2 
(4,)0 به صورت یک انتگرال محیطی قابل حل 


می‌باشد [6] حل عددی رابطه‌ی (۶۸) با استفاده از 


روش تجمع محلی به صورت 


1 1 (۱18,7 
(-- +2) ۳ 22 (4062--- سح 
۷ 7۳۳ 9 9 
بت ۳ 
2>۱ م((۷)(-)[(--- )و9 + 
۷ ۷۴ 


می‌باشد. می‌توان نشان داد تنش پرشی 2 روی 


مرز ۲-3 به صورت زیر محاسبه می‌شود: 
رم 0 (2 ,)وب 
- 9 
)2( ت 
۳ م ۳ ۷ 1 
[(ه + 12+ (2-6) ,16 + -] .۲ 
74 


۳( 


با استفاده از روش تجمع محلی برای حل انتگرال 
(۵۳) می‌توان رابطه‌ی اخیر را به صورت عددی 


نشریه مهندسی عمران دانشگاه فردوسی مشهد 


1۲ 


ِ رمرم 90 - (82)ورت 
(,۱۷,۲)۷ 2 +(2 »4 ۱ 


2-1 


(۵4 2 1 ِ ِ * 
۷+2۸ 2-1/2 ۷ 
ی 1 1 
ار )و( ه) )رت 
/ ۷ ۷۷ 
محاسبه نمود. 


برآورد عددی 


به منظور بررسی پاسخ محیط. تنش پیچشی متمرکز و 
معلوم (2) ۲ را در عمق 5 به صورت زیر تعریف 
می‌کنیم: 
(۵۵) 1 >0 ,(5- 2002 < (2) ۲ 
در رابطه‌ی فوق (۵)00 تابع دلتای دیراک می‌باشد. 
تغییرمکان ناشی از نیروی پیچشی 0 روی مرز 
۲<8 در عمق 920.5 برای مصالح همسان در شکل 
در 5< 2 بی‌نهایت می‌شود و با دور شدن از 5ع< 7 
به سرعت کاهش می‌یابد. در این شکل نتایج حاصل از 
این مطالعه با نتایج ارائه‌شده توسط یک و عابدزاده 
(۱۹۹۲) مقایسه شده است. تطابق عالی جواب‌ها معرف 


صحت نتایج به دست آمده می‌باشد. 


0 02 04- 06 08 
0 
2 
سس ت 
15 سا 4 
۱ 
1 
50 1715 
۳ (1992) 02۵۵602506 « 
15 ۱ 
7 ۳ 
--2 وا 
0 ۱ 
۱ 
ردو 


شکل ۲ تغییرمکان دیواره‌ی حفره در مصالح همسان ( 2 <1) 


برای بررسی اثر میزان ناهمسانی به رفتار محیط. 


علاوه بر محیط همسان. محیط‌های همسان جانبی با 


تشریه مهندسی عمران دانشگاه فردوسی مشهد 


تحلیل محیط نیمه بی نهایت همسان جانبی به پیچش وارد بر... 


ضرایب 0<0.5 و 1.5-<0 در نظر گرفته می‌شوند. 
برای محیط همسان نیز 1< »0 می‌باشد. شکل (۳) میزان 
جابه‌جایی دیواره‌ی حفره بر اثر اعمال تنش پیچشی 
(2) ۲ مطابق (۵0) روی مرز ۲<۵ و در عمق 21 ٩‏ 
را نشان می‌دهد. میزان تأثیر ناهمسانی مصالح بر 
جابه‌جایی به وجود آمده بر دیواره‌ی حفره به طور 


واضح مشاهده می‌شود. 


ی .مق 0.5 < 07 ,(ا) امعصععدام‌وزه - 

21 0 ,() ۳060۶عع2اموز0-- 
۵-2 215 ۵7 ,(د) احعطهعداموزه + 
2 3 


شکل ۳ تغییر مکان دیواره‌ی حفره مصالح همسان و 


همسان‌جانبی( 2 - 1) 


5 2 ۵ و(۲) 5۱۲۵55 5۱62۲ - 
2 0-1 ,(2) 5و5 وهی + 


215 ۵ ر(۲) 5655 5062۲ ۰ 


8 ,/ 
: 5 
شکل ۶ تنش برشی دیواره حفره مصالح همسان و 


همسان‌جانبی( 2 1) 


میزان تتش به وجود آمده در جداره‌ی حفره بر اثر 
اعمال تنش پیچشی (2) ۲روی مرز ۲-8 و در عمق 
21 و در شکل 63 مشاهد می‌شود. همان‌طور که 
ملاحظه می‌شود. تنش برشی به وجود آمده در دیواره‌ی 
حفره در محل اعمال تنش پیچشی برای همه‌ی 
مصالح. به سمت بی‌نهایت میل می‌کند. و با دور شدن 


از محل اعمال تنش پیچشی. مصالح رفتار متفاوتی از 


سال بیست و چهارم شماره یک »۱۳۹۱ 
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خود نشان می‌دهند و با نزدیک شدن به مرز 7-0 
میزان تنش برشی برای هر دو مصالح به سمت صفر 
میل می‌کند. میزان تنش برشی روی دیواره‌ی حفره 
برای مصالح متفاوت. در مرزهای 2-0 2-1 و در 
محل اعمال گشتاور پیچشی یکسان و در سایر قسمت- 
ها متفاوت می‌باشد. تفاوت تنش برشی در محیط‌های 
متفاوت مطابق شکل (4) به معنی اهمیت در نظر 
گرفتن اثر اایزوتروپی محیط در بررسی این رفتار 
می‌باشد. 


نتیجه گیری 
در اين مقاله یک روش تحلیلی برای به دست آوردن 
پاسخ محیط نیمه بی‌نهایت با رفتار همسان جانبی به 
پیچش وارد بر جدار حفره با طول محدود ارائه شده 
است. با تقسیم محیط به دو ناحیه‌ی ,3 و یل و 
اسفاته از ینبل کسوس فوزیه تمادلات تسادل زد 
دستگاه مختصات استوانه‌ای (۲,0,2) برای هر ناحیه 
به معادلات دیفرانسیل معمولی بر حسب عمق تبدیل 
شده به طوری که با حل آنها» شرایط پیوستگی دو 
اش ای رتیت وله منت 


مراجع 


1۳ 


حل معادله‌ی انتگرالی کوشی با استفاده از بسط تابع 
مجهول به چندجمله‌ای ژاکوبی و نیز استفاده از روش 
تجمع محلی» یک دستگاه معادله آ مجهولی 
حاصل گردیده که با حل آن‌ها تابع تنش برشی و سپس 
تابع تغییر مکان روی دیواره‌ی حفره در فاصله‌ی 
1> > 0 به صورت انتگرال‌های خطی نیمه‌متناهی به 
دست آمده است. توابع زیر علامت انتگرالی به علت 
وجود توابع بسل اصلاح شده توابع نوسانی می‌باشند. 
که در بی‌نهایت به سمت صفر میل می‌کنند؛ اما روند 
همگرایی این توابع در بی‌نهایت کند است. بنابراین؛ 
یکی از موارد مهم در برآورد انتگرال‌ها در این مقاله 
تعیین بی‌نهایت فیزیکی است. برای بررسی نتایج تنش 
برشی معلوم روی دیواره‌ی حفره در عمق ۶ فرض 
شده است و در نتیجه توابع تغییرمکان و تنش به 
صورت توابعی از 5 محاسبه گردیده‌اند. به منظور نشان 
دادن تأثیر میزان ناهمسانی مصالح بر پاسخ محیط نتایج 
مختلف برای تغییرمکان و تنش دیواره‌ی حفره ارائه 
گردیده اشتت: همان گونه که ملاحظه می‌شود به علت 
تفاوت‌های موجود بین نتایج مصالح همسان و همسان 
جانبی, لازم است در کاربردهای مهندسی برای نزدیک 
بودن پاسخ به واقعیت؛ مصالح به درستی معرفی گردند. 


و25 ۷۵۱۰ ,]۷6۵ .۸ .[ ,۲ ۱۵206 هطذمصط وه هام اقط ‏ صد 0ععنا0م۲ص فعوفع‌تاو. ر.ل رطاتا۲ 2۳0 . 0016 
1958(۰) ,443-436 .00 

,1-86 .00 ر4 ۷۵۱ ,1000 ۷6۵۰( ," فتصهطعهه صا فصماطميص مضه ۱۷۲۵0 ربط رطقم۳۲00 
.(1979 

1944(۰) ر00۲طم۱ رقفعزص اتفتع تا عمتصصصه ر< بچا‌تادهاع ۰۲۳۱6۵۲6۵6۵1 ,۷۷۰ ر26۳۳98 200 ۳۰ .ظ رطعع۲) 
یک ,۷۵ ,۷6611 ۷0 .ل ,۲ 5۵۱۲۵ ۷۲6 نوماه آم‌تسنان مانصت ه مت ۷۷۵۲۵ فوعتاو ۳۲6 ,12۰۲۷۷۰ ر0۲02۳[ 
1962(۰) ,503-551 2۳۰ 

(1981) ,۱۷0۹60۲ رتعطوناصانظ نا ر" رقم »زوم هدنج 6 ]۵ وت‌تاعهاه ]0 ۰1۳6۵ ری.9 رتلاکاتصطاع,1 

و 21۲ امعنصنان عاتصت ۵و۵ جح صم ماع اقعمن108 ر.۲ رطاع02۵0عظ۸ 220 ,۴۱۷۰9 رعا۳2 
.(1992) ,133-144 .00 ر438 ,۷۵۱ رش .10۳0 506 

رکه 5۳۱۱60 5015 . ۲۲ ,۲ 90۲6 آ‌نهصنان عاتصتص نو ۵۴ ع۲۵6] و عظ) وم فصمتا۵ ۰۸۵۵۲60 ربظ ر5ع2۲۳۲6ظ 


7۵۱, 19, 2۲۰ 165-177, )1982(, 


سال بیست و چهارم شماره یک ۱۳۹۱ نشریه مهندسی عمران دانشگاه فردوسی مشهد 


34 تحلیل محیط نیمه بی نهایت همسان جانبی به پیچش وارد بر 


0 و" 268و 0۵۲6 ۵ جه فصتاعه ع0۲6] موم »تطمصتقطعصتا و ما عفع0حوع1۲ عتافقاط رب ر5ع ۳210‏ .8 
1983(۰) ,925-934 00۵۰ ,19 ۷۵۱۰ رکع 506 علز501 

جر و ۲اه یسنان و ۵۶ عصندم۱ عساووعتج 0عاحتاجععصم ۵ 1۵ متا فمز)زتمانامطله ۵ بظ رفعصتوظ .. ,9 
.(1984) ,267-276 00۰ ,20 ۷۵۱۰ ,50۳695 عل[501 

و 5۲۱610۳65 50105 . ۲۱ ,۲ تلع تاعقاع جع صز عت۵ه ‏ هصماه و۱۵۵ [مصمنعتما مصلقومتم۰۳۳۵ رب رقع2۳0ظ .10 
.(1980) ,653-670 .00 ر16 .۷0۵1 

۰ (1۲6۹ع 5۳ 50105 .ل 1 ,۲ 0266و عمط 0۲۵ و ج۵ عتاووعتم ۲1۵-1020 عاهاو-۱6۵0۲" رب ,۳27۳68 ,11 
1986(۰) ,73-86 00۰ ,22 

و ۳۷5۰ ۸۳۷ ۲ فعازی ام‌زهستان 4و ای‌نمطاوی مر ۲۵۲۵۵ ممتوفع۲مصم ]۲2۳۲6160 ]۰ ربا ۲۷۷۰ ,5610678 .12 
<(1952) ,97-108 0۳0۰ و5 .۷0۵1 

,(1951) رطمتان۴ 1٩‏ رتم۷ مه :۱۷60-۳۳111 رد عوم کم م سم ربلا:1 رح080عع5 ,13 

عاتگن هو ۵ ۵۶ ۲26۵ عطا عصمله هتم عساععهتج ۵۶ عولنام ج 0۲ ۳۵01۲680( وفع‌تا" ,.ل:1 600ص 14 
.(1952) ,57-62 0۳۰ و2 ۲۵16۲0 نل ممتامم۱۷۲6 مامت 6۵060 ,50110 

1941(۰) ردل)عالهن ر134-161 .20 رعصنا0 ۷ 7۵۳5۵۲۲ تحص مقصصتق! مب ۰۲۳۵۵00۲۵ ربک۲۲:1۷ ر۵۵۲6۲۵88۲ ۷۷ .15 

1. خجسته ع. رحیمیان. م و اسکندری قادی. م. «تحلیل سه‌بعدی محیط نیمه بی‌نهایت با رفتار ایزوتروپ جانبی تحت اثر بار مماس بر 

سطح در فضای فرکانسی». نشریه‌ی دانشکده‌ی فنی» جلد ۶۰ شماره ۵ ابان (۱۳۸۵). 
۷ فرشاد. م. «ریاضیات عالی مهندسی»»جلد دوم انتشارات بعثت» چاپ اول» (۱۳۱۷). 
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